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Abstract 

We present the perturbative matching coefficient to 0(aa s ) which relates 
the AB=2 operator in the continuum to that of the lattice static theory, 
which is important in the accurate extraction of the continuum value of the 
Bb from lattice simulations. The coefficients are obtained by the one- loop 
calculations in both of the continuum and lattice theory. We find that two new 
dimension seven operators appear at the 0(aa s ) with the 0(1) coefficients. 
We also discuss the possible cancellation of 0(aa s ) correction in the ratio 
B B = (B\0 L \B)/((8/3)(f B M B ) 2 ) qualitatively. 

PACS number(s): 12.38.Gc, 12.39.Hg, 13.20.He, 14.40.Nd 



1 



Typeset using REVTgX 



I. INTRODUCTION 



One of the most important issues in particle physics is the origin of masses and CP 
violation. CKM matrix elements are believed to play a key role to probe the physics behind 
it. Despite a lot of efforts in various approaches, the matrix element V t d which can be 
determined from B° — B° mixing is still only poorly known because of a lack of accuracy in 
the involved hadronic matrix element. The hadronic matrix element is parameterized using 
the B meson decay constant fs and the bag parameter Bb as Bb]\i so ^ * s quite crucial to 
compute them with high precision. For this purpose, the lattice QCD has been considered to 
be one of the most reliable approaches. So far most of the efforts have been devoted to the B 
meson decay constant. At the early stage, the decay constants in the static approximation 
and from the extrapolation from light quarks were computed. It was found that both the 
lattice cutoff dependence and heavy quark mass dependence are significantly large. Later 
the scaling behavior for the lattice spacing a |TJ|^] and the heavy quark mass 1 / irtq || have 
been investigated carefully and the best estimate on f B from the quenched Lattice QCD is 
now Jb = 165(20) MeV Q. On the other hand, until recently, the bag parameter has been 
calculated only either in the static limit or by the naive extrapolation from light quarks. In 
this respect, careful studies of systematic errors of the bag parameter are still missing. 

In general, in order to get a continuum result of the physical quantity such as f B from 
lattice simulation, we have to compute physical quantities on different lattices and extrapo- 
late the results to the continuum. Therefore the final results have smaller errors if the cutoff 
dependence is smaller. It was found that the 0(a) improvements of the action and lattice 
operators in Symanzik approach significantly reduce the lattice cutoff dependences of various 
matrix elements. For the heavy-light axial vector current, such kind of improvements have 
been accomplished by Morningstar and Shigemitsu in the lattice NRQCD formalism. 
They found that the additional operator mixed at the 0(aa s ) and the inclusion of the effect 
significantly reduced the value of fs at the finite lattice spacing and it was also the case in 
the static limit. In contrast to the decay constant, the 0(aa s ) mixing effect has not been 
studied for Bb- One reason is that only the operator matching of 0(a s ) has been done in 
Refs. |6|-§3| so far. Although previous simulations have not shown a clear cutoff dependence 
of the B B [@,[HJ, it would be very important to study the 0(aa s ) mixing effect explicitly in 



order to obtain the precise value of Bb- 

The purpose of this paper is to investigate the 0(aa s ) effect for the Bb- We perturba- 
tively compute the operator matching coefficients of the static-clover AB=2 operators up 
to 0(aa s ). We use the notation defined by the authors in Refs. ||[7||. 

Phenomenologically important quantity might be the product of BbJb which is just 
the expectation value of AB=2 operator. Therefore it seems sufficient to improve only 
the AB=2 operator. To determine /# and Bb separately, however, would have some more 
advantage from a technical point of view [|]]. Since the 0(aa s ) improvement for Bb requires 
the improvements of both of the heavy-light axial vector current and the AB=2 operator, 
we also mention the result for the heavy-light current for completeness. 

The paper is organized as follows. In sections [IT] and |T|, our main results, the matching 
coefficients to the 0(aa s ) for the heavy-light current and the AB=2 operator, are shown, 
respectively. In section [TV], we discuss the impact of our results on the determination of 
the B B using the typical values of the parameters involved and the consistency with the 
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previous observations for the cutoff dependence of the Bb- Finally we conclude in section 
[V[ The appendices are devoted to some details in this calculation. 

Throughout this paper, we choose Feynman gauge (a=l) and the light quark mass m q 
is set to be zero. The ultraviolet divergences appearing in the continuum calculation are 
regulated by dimensional regularization and the continuum operators are renormalized with 
MS scheme, while the infrared divergences are regulated by the gluon mass A in both of the 
continuum and lattice theory. The operators with superscript "con" and "lat" define the 
continuum operators and the lattice operators, respectively. In our convention, 75 always 
anticommutes with 7^. We give all the equations in Euclidean form. 



II. STATIC HEAVY-LIGHT CURRENT 

In this section, we present the matching coefficients of the static-light current operators 
which are relevant to the determinations of the form factors of the static to light decays as 
well as the following discussion. Our lattice gauge action is the standard Wilson plaquette 
action. For the light quark we use the 0(a)-imp roved SW quark action |TT| with the clover 



coefficient c sw and, in contrast to Ref. J/J, we do not incorporate the rotation operator 
associated with the clover fermion in the current operator. 

In the following, we describe the lattice static quark. In the static limit, the quark action 
is separated into two pieces in Dirac basis, namely one for the static quark b' and the other 
for the static antiquark b'. Both are two-component fields which are related to the relativistic 
four- component field b as 




(&* 



(1) 



In our convention, the action is given by 

+ E (-vUx)) I V - t?WWM?(»), ( 2 ) 



x,y - ' 

where a (a') and (3 (J3') run over 1 and 2 (3 and 4). Our Feynman rules for the lattice static 
quark and antiquark are obtained from the above action through the standard procedure. 
The heavy quark (antiquark) propagates only forward (backward) in time direction on the 
lattice. 

To determine the matching coefficients upto 0(aa s ), (i) we calculate the heavy to light 
on-shell scattering amplitudes through the following operator with arbitrary gamma matrix 

r, 

in the continuum full theory upto the one-loop order, expand the resulting expression with 
respect to the momenta of external quarks at their rest frame, which is required to obtain 
the matching coefficients through desired order 0(aa s ), and take the static limit of the 
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heavy quarks, (ii) We repeat the similar calculation to step (i) on the lattice static the- 
ory, (hi) Finally we express the continuum operators in terms of the lattice operators with 
appropriate matching coefficients which are adjusted to coincide the one-loop scattering am- 
plitude of both theories through 0(aa s ). In this matching procedure, we have two coupling 
constants, in the continuum theory and aj, at on the lattice theory. Through this paper, 



T3J at one-loop 



both coupling constants are rewritten in terms of the V^-scheme coupling 
order. 

According to step (i), we calculate the scattering amplitude with an initial heavy quark 
carrying momentum p and a final light quark carrying momentum k. The resulting expres- 
sion is 



( q (k)\ 4 0)con Hp)) 




(3) 



where the symbol (• • -)o denotes the tree level expectation value between the same initial 
and final states as those of the left hand side, Cf = (N% — 1)/2N C with number of color 
N c , mb is the heavy quark mass, and Jp = q(aD ■ j)Tb. The renormalization scale for 
the amplitude is /i. The definitions of H, G and H' are the same as those in Ref. |14" |. In 
deriving Eq. we use the equation of motion for the light quark, qj^k^ = —qj ■ k and 
also that for the heavy quark, 7 4 -u 6 = w b , to simplify the result. 

Repeating the similar calculation to the continuum theory according to step (ii), we 
obtain the corresponding amplitude on the lattice as follows. 



(g(*)| 4 0)lat Hp)) 



1 + ^° F (-^n(« 2 A 2 ) + 4 0) +4 (0) + ^ 2) 



/ t(0)\ 

Wr /o 



G 



in 



3a\ 



+ r (1 - c s 



^AVAf'+if^jfto, (4) 



where 



.4 



(o) 



.4 



1(0) 



A 



(i) 



.4 



di + d 2 G + -[e 

-d'G+^f 1 , 

UG + V, 
U ! G + V 1 . 



+ f), 



(5) 

(6) 

(7) 
(8) 



The renormalization scale for the amplitude is a 1 . A^ and A^ correspond to Ar and 
A 1 in Ref. J7|, respectively, and the numerical values of di, d 2 , and / are tabulated in 



Refs. P,[7jn|. Although our explicit form of the integrand of d 1 completely agrees with that 
of Ref. Jj]], the numerical value of d 1 is slightly larger in magnitude than that of Ref. 0, 
and the value is tabulated in Table |. U, U 1 , V and V 1 are new contributions at 0(aa s ). 
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Their explicit forms of the integrands are shown in Appendix [B] and their numerical values 
are tabulated in Table Q. The coefficients with the superscript I vanish when Wilson light 
quark is used (c sw =0), which is the same notation as Ref. [0. comes from the tadpole 
improvement of the light quark wave function renormalization, for details see Appendix [A|. 

In step (iii), matching Eq. (|3]) to Eq. (^), we obtain the following relation between the 
operators in the continuum and lattice theory, 



(O)con 



i + ^c F c r 0) 



Air 

z (0) j(0)lat + z (l) { C }] 



7 (0)lat 0! s (i) (i)i a t 

j r -+- sr J r 



Sr U 2 J [mi 2 [a 2 2 4 4 



where 



+l\n(a 2 \ 2 )-A\"' -.4;."" (10) 
C W = _ r (1 _ Csw) ln(a 2 A 2 } _ A », _ A p, (n) 

Eq. (|10|) has been calculated in Refs. 0,0] except for the differences of our inclusion of tadpole 
improvements and the wave function renormalization of lattice static quarks. Eq. ( |Tl~D is new 
result for the arbitrary static-light current. For axial vector current and vector current the 
matching coefficient for Jp^ lat has been calculated with NRQCD action for heavy quarks in 
Ref. [Q. From Eq. (|^) we observe that the 0(a) operator Jp 1 ' )lat appears at this order, which is 
considered to be a lattice artifact. It is noted that there is no linear divergence proportional 
to 1/A in the coefficients, while there is a logarithmic divergence unless c sw = 1. In the use 
of Wilson light quark (c sw =0), therefore, we cannot match these operators consistently due 
to this infrared mismatch as previously pointed out in Refs. [0,0. 

The results of and for each T are summarized in Table [0], where r = c sw = 1 and 
the tadpole improvement is performed by using the perturbative expression of the critical 
hopping parameter. The numerical values of 0(aa s ) correction for axial vector current and 
vector current are consistent with those in Ref. ||[]. It should be noted that the coefficient 
of J r 1 ^ lat depends only on G, G = — 1 might lead to a large mixing effect, while G = 1 does 
not. Actually the mixing effect leads to the significant change for the fg, which has been 
seen in Refs. [@,^]- 

III. AB=2 OPERATOR 

In this section, we discuss the matching of the AB=2 operator. The matching procedure 
is completely common as that for heavy-light current previously shown, and we follow the 
previous step. Before proceeding step (i), we give the definitions of the operators. 



Note that since there are some differences of the definitions of the lattice operators and the 
matching coefficients between ours and theirs in Ref. ||, one would need to redefine our definitions 
to compare the results with theirs. 
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O s 

Or 

O n 
O ld 
Ond 



bP L q] [bP L q] , 

■ \fhi*PLq] \h»PRq\ + 4 
h^PLq] \h t iPR(aP ) ■ i)q 



bP L q\ [bP R q 

+ A\bP L q bP R {aD-^)q 



where Pl=1 — 75 and Pr=1 + 75. 

According to step (i), we calculate the two-body scattering amplitude for Ol between 
the initial state with a heavy antiquark and a light quark and the final state with a heavy 
quark and a light antiquark in the continuum theory. The initial heavy antiquark carries 
momentum ~p 2l the initial light quark k 2 , the final heavy quark pi, and the final light 
antiquark k\. We obtain the scattering amplitude in the continuum theory at one- loop as 

(qfyxM or \qfyxp2)} = z™zrY,v$Api,k2,P2) 



a 



mt 



l + 2 In -£ -4 In ^ \+C L + - 



Att 



mt 



(Ol)o 



+^C s (O s ) + ^^(O ND ) , 
An Att SaX 



(12) 



where the V$ n (i runs over a-d) denotes the contribution from each diagram in the continuum 
theory, which appear in Appendix |C[ The constants Cl = — 14 and C$ = —8 appear in 
Refs. |,| . 

According to step (ii), we calculate the corresponding amplitude with the lattice theory 
and obtain the result as follows. 

(qfyMP^l OT \qfy,b(P2}) = Z^Z^Y,^(ki,Puh,p 2 ) 



1 + g (-41n(a 2 A 2 )-D L -D[ + ^ + ^) 



(Ol)o 



+^ {-Dn ~ D N )(O N ) + g (-D R - D R )(O R ) 



+ ^ ( r(l - c sw ) ln(a 2 A 2 ) - D LD - D[ D ) (O ld ) 



+ 



An 
o.s- 
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where 



3 

/16vr 
4^ \SaX 



D, 



Dnd — Dnd ) (Ond)o, 



10 , 1 1 

Cti c V 

3 3 3 



+ f) + 



D 



N 



V 

3 

2d 2 , 



(13) 

(14) 

(15) 
(16) 
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Dn = 


-2 d 1 , 


(17) 


n 


4 

-to, 
3 




D R = 


4 , 


(19) 


D L d = 


3 


(20) 


Dld = 




(21) 


Dnd = 


-2U, 


(22) 


D 1 - 

U ND — 


-2U 1 . 


(23) 



The coefficients D L , D L , D N , D ! N , D R , and have been calculated in Refs. and we 
use the same notation as those in Refs. for convenience. The coefficients Dld, D ld , 
Dnd, and D^p are novel results of this paper. V^t (i runs over a-d) is the contribution 
from each diagram in the lattice theory, which are shown in the Appendix |C[ 

According to step (iii), using Eqs. ([12]) and flT3|) we match the the lattice operator and 
continuum one to 0(aa s ). The obtained operator identity is 

O c r = J2ZxOf, (24) 

x 



where X runs over {L, S, N, R, LD, ND} 
Z L 



1 + ^6 \n(a 2 m 2 b ) - 2 ln(a 2 /i 2 ) + Cl + D l + D[- -u^J , (25) 

Z s = y- C 8 , (26) 
Z N = ^ (D N + D T N ) , (27) 
Z R = ^ (Z>* + , (28) 

Zld = £ ( y r (1 " Csw) ln(fl2A2) + Dld + ' (29) 

Z WD = |i (l^ D + Z^ D ) . (30) 

Here we omit the explicit arguments of /x and a -1 , which are introduced by the renormal- 
ization procedure, for the operators Ox and the coefficients Zx without ambiguity. We find 
that the above results to 0(a s ) agrees with those of Refs. except for the coefficient 
in Ref. M (see Appendix y). The correct value of including double rotation operator 
has been already obtained in Refs. |||| and our is consistent with them. Two new oper- 
ators Of D and Offiu mix at the 0(aa s ). It should be noted that the coefficients of the new 
operators have completely common integrands to those of Jp 1 ' )lat in the heavy-light current. 
The use of Wilson light quark (c sw =0) leads to the mismatch of infrared behavior between 
continuum and lattice theory as in the case of heavy-light current. 

When c sw = r = 1 and the tadpole improvement by the critical hopping parameter are 
chosen for the numerical estimate of Eqs. (E5])-(B0), Eq. (E4j) becomes 
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o 



con 



l + ^-U Ma 2 m 2 b ) - 2 ln(aV) - 35.15 



O 



lat 



+^ (-8) 0^ + ^ (-6.16) 0* + -± (-0.52) 0* 
+^ (-17.20) Ojg, + ^ (-9.20) 



(31) 



It is found that the coefficients of two new operators are 17.20/47T and 9.20/47T, respectively, 
and are of 0(1). This means the possibility of large 0(aa s ) correction for Of n as in the 
case of axial vector current, though the lattice matrix elements of Of D and Of D are not 
yet known. 



IV. DISCUSSION 



In the previous section, we pointed out that the AB = 2 operator might receive a large 
0(aa s ) correction. For the rigorous investigation of the 0(aa s ) effect, we must rely on the 
future works. On the other hand, the previous simulations have not shown a clear cutoff 
dependence of the B B and seem to imply that the vacuum saturate approximation (VSA) 
is plausible within 10% level around the used lattice cutoff scale (~ 2-3 GeV) P, |lC||i6| . |T7| . 
In this section, therefore, we attempt assuming the VSA for the lattice matrix elements to 
estimate the 0(aa s ) effects for the B B f B and B B using the results of the previous sections 
and then investigate the consistency of our result with the previous simulations. Although 
this analysis is quite rough, we believe that it is possible to find some, at least, qualitative 
features. 

Let us discuss the 0(aa s ) correction for (B°\ 0™ n \B°) using the VSA. Under the VSA, 
the relevant lattice matrix elements take the following values, 



(VSA) 



(B°\ Of |5°) {VSA) = (B°\ Of \B°) 
(W\ Of |£?°)( VSA ) = --(W\ Of |5°)( VSA ) = - (f 



c (0)lat 
' B 



Mi 



(B°\ Of D |£°) (VSA) = (B°\ Of D \B°) 



(VSA) 



(32) 
(33) 



where f { B )Ut M B = (0\J^f\B°) and 8f\ 
Eqs. (|32| ) and (|33"D into Eq. ([H]), we obtain 

(W\ 0™ n |j?°) VSA > (B~o\ O c ^ n |S°)( VSA ) 



lat 
IB 



(0 1 jW'f \B°)/(0\ mf \B°). Substituting 



l + £ 61n( 



2 2\ 

a m h ) 



2 ln(a V) - 36.83 + 26.40 5f% 



lat 



(34) 



where the last term with 5ff is essentially due to the 0(aa s ) effect. We can use the data 



calculated by Ali Khan et al. in Ref. ||T5[ to guess the value of 5ff in the static limit. In our 
estimate, their finite mass results at (3 = 6.0 imply 5ff ~ —0.5 in the static limit. Using the 
coupling constants at the corresponding lattice with Lepage and Mackenzie prescription [ I3 |, 
a s ~ 0.15-0.25, we find that the magnitude of the 0(aa s ) correction for (B°\ O™ 11 |£? )( VSA ) 
is very large, about 15-25%. Although this analysis is naive estimate of 0(aa s ) correction 



S 



using the VSA, this suggests that there are large contribution from 0(aa s ) correction for 
Of 11 and the improvement of 0(aa s ) should be necessarily included. 
Now we turn to the Bb- The Bb is defined by 



B 



B 



(b°\ e>2 on I b°) 



UbM, 



(35) 



B 



To improve the Bb in consistent way, we should include the 0(aa s ) improvements of both 
of the numerator and denominator of Eq. ([ffi). Substituting Eqs. (JD and fl24|) into Eq. (|35"D 
and linearizing the resulting expression in a s according to the discussion of Ref. [|1(J, we 
obtain the Bb as follows. 



x 



where X runs over {L, S, N, R, LD, ND}, 



U X 



( 7 (0) \2' 

7574 
-(0) ' 



^7574 



B 



lat 



(B°\ Of \B°) 



In the VSA, using Eqs. (B2H and (Rl) we obtain the following expression for B B to 0(aa s ). 



BbY sa fl cvsA) 



' B 



u l + u r + uj n - -u s J - (uj ld + u ND + 2 u x ) 5f 



lat 
B ) 



(36) 



where the term with Sfg comes from the 0(aa s ) improvements again. The coefficients D 
and E are given as follows. 



D = [21b ^ - ^ - I (dx + d 2 - d 1 ) - \c- ^v + v 1 ) + t (w + w 1 ) 
r (1 - c sw ) ln(a 2 A 2 ) + C/ + U 1 + V + V 1 . 



2 
3 



(37) 
(38) 



In deriving Eqs. (|37|) and (p8|), there are some cancellations between the coefficients of the 
AB = 2 operator and the axial vector current. 

Now let us roughly estimate the 0(aa s ) effect in the B^^ numerically. When r = 
c sw = 1 is chosen, we obtain D = 21n(m 2 y / u 2 ) — 3.72 and E = —0.37. Using the data of 



Sf^ and the coupling constants as before, we find that the 0(aa s 



effect for the B^^ is 

smaller than 1%. Of course such a drastic cancellation would not take place in reality due to 
deviations from VSA, but at least the present analysis suggests that there is a possibility of 
significant cancellation of 0(aa s ) corrections in Bb- This is consistent with the observation 
from the previous simulations that there is not clear cutoff dependence in Bb- 
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V. CONCLUSION 



In this paper, we reported the coefficients of the 0(a) operators which are newly induced 
at the 0(aa s ) in the perturbative continuum-lattice operator matching of the heavy-light 
current and AB = 2 operator. We also roughly estimated the 0(aa s ) effect on the B B f B and 
the Bb using the VSA in the lattice hadronic matrix elements. Although the 0(aa s ) effect 
is significant in the determinations of /# and BbJbi it seems that the effect is not so for 
Bb, at least, in this VSA analysis because the cancellation between the 0(aa s ) effects in the 
numerator and denominator does work well. Therefore the previous works, which imply that 
there is no cutoff dependence in Bb, seem to be consistent with our analysis. Now, however, 
that the 0(aa s ) improvement for the fs has been already done, in order to calculate the 
Bb in a consistent way the 0(aa s ) operators should be included in the calculation. For 
the precise determination, it is also required to include the finite mass correction into both 
calculations of the matrix element and the matching coefficients. 
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APPENDIX A: 

Here we show the wave function renormalization constants for each external quark line 
in each theory, 



z con = j 



Z c b on = 1 



7 lat 
^ lat 




-21b 



A 2 



4n 



ln(a 2 A 2 ) + / + f + u (2) 



Zf x = Z, 



lat 



Hat 



4:71 




2 \2 



+ 4 



-21n(a 2 A 2 ) + e 



where /, f 1 and can be found in Refs. ||[7|. In above equations A = 1/e + ln(47r) — 7^ 

(2) 

and Uq is perturbative coefficient of the tadpole improvement factor defined by uq = 1 + 

(2) (2) 

ch s CfUo . The coefficient u is obtained through the calculation of the mean plaquette value 



or the critical hopping parameter, 
respectively. 
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(2) 



-7r or u 



(2) 




[4.4259 + 8.4327 r - 4.8619 c s 
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APPENDIX B: 



Here we show the explicit forms of the integrands for U, U 1 ,V and V 1 , which first appear 
in the 0(aa s ). For shorthand notation, we define the following quantities, 

4 



A 1 = 



Ao = 



Af = 



^sin^g+^^AO 2 , 
/i=i 



Af = £sin 2 (y+4r 2 (Af) 2 , 
/i=i 

3 



3 // ' 

2 / l /u 



^sin 2 (/^) sin 



Using the above convention, 
r* d 3 l 



12Af Af 



(3 + (3r 2 - l)Af 



(3) 



1 



(Af - 2Af + 2r 2 A^Af ) - 



(3) 



2 A (3) A (^ 



12Af (Af ) 2 



-0(1 -P) 



3(Z 2 ) 



16 

y 



^ = (4vr) 2 r 2 f 

J — IT 

V = (4vr) 2 r 



ci 3 Z 
(2tt) 3 

d 4 / 



.(3) 



(Af - 2Af + 2r 2 A^A^' 



.(3) 



48Af Af 12(Af) 2 
1 



2 A (3)a(3)> 



4A 9 



- 12Ai X (A2)2 {12 (l + 2Af + 2(r 2 - l)Ax) (1 - A 1 + Af )A?> 
+(Af - 2Af + 2r 2 Af AO} + ^0(1 - / 2 ) 



(l 2 ) s 



/■7T ^4; 



(2tt)< 



12Ax(A 2 )2 



{ (l + 2Af + 2(r 2 - 1) Ai) A 
+(Af - 2Af + 2r 2 Af Ai)} (1 - A x + Af ) - ifl(l - Z 



(3) 
4 



APPENDIX C: 

Here we show the contribution from each diagram explicitly. In the continuum, the each 
contribution is as follows. 
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y(«0 

r con 
K con 

r con 

v rc\r\ 



(Ol)o, 




11 

y 

2 

" a 

5 
3 



(Ol)o 



4^ 8 ^ + 4.3aA 



H 167F /ro \ 



In 



mi 



(Ol)i 



(Ol) 



0- 



where (a) corresponds the tree diagram, (b) those with the gluon connecting the static and 
the light quarks, (c) those connecting the static quark and the static antiquark, and (d) 
those connecting the light quark and the light antiquark. And on the lattice, 



V, 



(a) 
lat 



K lat 



(O L ) , 
a, 10 



4tt 3 

a s 

+7^ 2 

10 r 



+ 



V (c) 

Hat 
K lat 



\~ix ~3 
+ 7^ 2 

47T 

1 

4tt 3 

Og_ 1 

4tt 3 

« s 4 
+ — - 
4tt 3 



-huVA^ + dxj (O L ) 
-da + d 7 ] (Ojv)o 

r(c sw - 1) ln(a 2 A 2 ) — (V + V 1 )] (O LD ) 



87T 

3aA 



+ (U + U 1 ) 



21n(a 2 A 2 ) + c 



(Ond)o, 
(Ol)o, 



AHa'X^ + iv + v 1 )} (O L ) 
-(w + w 1 )} (O R ) . 



(CI) 



The calculation is straightforward though slightly lengthy. The full use of the equations of 
motion for the heavy and the light quarks and of the identities for 7 matrices sometimes 
leads to simplification, in particular for the derivation of V^ . We find that our result of 
V^£ is inconsistent with Eqs. (B.16) and (B.25) of Ref. [[7] provided the sign of the numerical 
values tabulated in Table 3 of the reference was correct, which has been already pointed 
out in Refs. MM. 
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TABLES 



TABLE I. The numerical values of d 1 , U, U 1 , V and V 1 for each value of r. 



r 1.00 0.75 0.50 0.25 

d 1 -4.14 -3.74 -3.12 -2.04 

U 4.89 5.27 6.16 8.26 12.72 

U 1 -0.29 -0.11 0.02 0.06 

V -7.14 -7.51 -7.72 -6.99 

V I 1.98 1.82 1.51 0.98 



TABLE II. The results of the heavy-light current matching and H, H' and G for each L. 



L 


H 


H' 


G 


Cr° } 






1 


4 


1 


1 


3 fhV/mD + |ln(a 2 m2) 
3 |ln(/x 2 /^fc) + |ln(a 2 mf) 


- 2.25 


0.56 


75 


-4 


-1 


-1 


- 8.41 


9.76 


li 


-2 


-1 


-1 


|ln(a 2 m^) - 14.41 




9.76 


74 


-2 


-1 


1 


|ln(a 2 m 2 ) - 6.25 




0.56 


757i 


2 


1 


1 


|ln(a 2 m 2 ) - 8.25 




0.56 


7574 


2 


1 


-1 


|ln(a 2 mg) - 12.41 




9.76 







1 


-1 


- |ln(/x 2 /m^) + |ln(a 2 m^) 

- |ln(^, 2 /mg) + |ln(a 2 mf) 


- 14.41 


9.76 


<Tij 





1 


1 


- 8.25 


0.56 
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